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Abstract. In superconductors with large values of the Ginzburg-Landau 
parameter k, exposed to magnetic fields close to the upper critical field -ffc2i 
the magnetic field is practically homogeneous across the sample and the density 
of supercurrents is negligibly small. In this case, there is no obvious reason for 
the formation of Abrikosov vortices, characteristic for the commonly known mixed 
state. We consider an alternative model for describing the mixed state for k ^ 1 
and magnetic fields close to Hc2 ■ We argue that with decreasing magnetic field the 
traditional vortex structure is adopted via a first order phase transition, revealed 
by discontinuities in the magnetization as well as the resistivity. 
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It is commonly accepted that the mixed state of a type-II superconductor is 
characterized by the penetration of an external magnetic field into the sample along 
quantized vortex lines or vortices. The spatial extent of the superconducting order 
parameter ipii^) near a vortex axis is determined by circular currents flowing around 
the vortex line. The velocity of the relevant charge carriers, i.e., the Cooper pairs, 
diverges at the center of the vortex line (see, e.g., Ref. thus reducing the 

superconducting order parameter progressively until it vanishes on the vortex axis. 
In most cases the vortices may thus be considered as thin normal filaments embedded 
in a superconducting environment. A completely different situation may, however, 
be established in superconductors with a Ginzburg-Landau (G-L) parameter k ^ 1. 
If the applied magnetic field is close to the upper critical field Hc2, the distance 
between adjacent vortex cores is much smaller than the magnetic field penetration 
depth A(T). In this case, there is practically no expulsion of the magnetic field 
from superconducting regions and the density of shielding currents is negligibly small. 
These circumstances are naturally unfavorable for the formation of common Abrikosov 
vortices and hence the character of the mixed state in magnetic fields close to Hc2 may 
be very different from that adopted at lower fields. 

In this paper we consider that in magnetic fields close to Hc2, the natural 
alternative to the conventional vortex structure is the formation of superconducting 
filaments, embedded in the matrix of the normal metal. The properties of such 
superconducting filaments may be analyzed by numerically solving the G-L equations 

First, we consider a single cylindrical superconducting filament in an infinitely 
extended normal metal. The magnetic field inside the normal metal is oriented parallel 
to the filament and its value is set to Hq < Hc2- For the numerical analysis we have 
chosen cylindrical coordinates (r, 0, z) with the z-axis parallel to the filament and 
r = at its center. The G-L equations for an infinitely long cylindrical filament 

and 

d /dA _ ^2) 

dr \dr r J 

may thus be solved in one-dimension, with all quantities depending only on the radial 
coordinate r. Here, tp is the superconducting order parameter normalized by its 
equilibrium value in a bulk superconductor; A is the vector potential of the magnetic 
field expressed in units of \/2Hc\(T), where He is the thermodynamic critical field. 
The coordinate r is measured in units of A(r). 

The boundary conditions are set at the center of the filament (r — 0). For the 
superconducting order parameter i^{r), we have d^ / dr\r={) = and i/'(0) — const. 
For the vector potential A(r), we have A{Q) — and curl A|r=o — H{0), the value of 
the magnetic field at the center of the filament. We have checked that the solutions 
with A{0) — 0, equivalent to a vanishing supercurrent at r = 0, correspond to a 
minimum of the free energy. By integration of the G-L equations with these boundary 
conditions, we obtain ip{f)i ^(^) f^nd the distribution of the magnetic field H{r) around 
the filament. We also calculated the Gibbs free energy Ffu of the filament per unit 
length. For a single filament, Ffn reaches its minimum if, at the boundary with the 
normal metal, the order parameter vanishes with a zero derivative, i.e., ipii^) ^0 ^^'^ 
d^j/dr for r — s- oo. 
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Figure 1. Radial variation of (a) tlie normalized order parameter ^^{r) and (b) 
the magnetic field for different values of Ho/Hc2- = \{T)/k is the G-L 

coherence length. 



The calculations have been made for 3 values of k = 10, 30, and 100. Fig. 1(a) 
shows the calculated profiles of the normalized superconducting order parameter 'ip{r) 
for one filament and, in Fig. 1(b), corresponding profiles of the magnetic field have 
been plotted. For all k values between 10 and 100 and Hq/Hc2 > 0.3, these profiles 
are independent of k. In this approach, all characteristics of the filament are uniquely 
determined by the value of the applied magnetic field Hq. 
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Figure 2. Different parameters of a single superconducting filament as functions 
of (1 — Hq/Hc2)- (a) The amplitude of the normalized superconducting order 
parameter at the center of the filament, (b) The normalized difference between 
the magnetic field at the center of the filament and Hq , multiplied by . (c) The 
Gibbs free energy of the filament Ffn (per unit of length) in units of ^■^ {T)H^ /8tt 
with respect to that of the normal metal in a magnetic field equal to Hq. 

In Figs. 2(a-c) we display various quantities, plotted versus (1 — Ho/Hc2)- Figure 
2(a) shows the amplitude of the normalized G-L order parameter at the center of the 
filament; in Fig. 2(b) we have plotted the product k^[1 — H{0)/Ho], and in Fig. 2(c) 
we show the free energy of the filament Ffu. It may be seen that while ip{0) and 
Ffii are independent of k, the expulsion of the magnetic field from the center of the 
filament, [1 — H{0)/Ho], is inversely proportional to k^, i.e., the magnetic moment of 
the filament rapidly decreases with increasing k. According to Fig. 2(c) the free energy 
of the filament with respect to that of the normal metal is negative and it vanishes at 
Hq = Hc2- This obvious but gratifying result demonstrates that the superconducting 
phase can only be stable in magnetic fields below Hc2- 
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So far, we have considered a single superconducting filament in a normal metal 
matrix and we now consider the interaction between the filaments. Because the current 
density in the normal metal between the filaments is zero there is no long-range 
interaction between them and the contribution to the free energy due to filament- 
filament interaction is zero unless the neighboring filaments are very close to each 
other. In order to show that the short-range interaction between the filaments is 
repulsive, we introduce the velocity of superconducting electrons Vg — js/{2eip'^), 
where js is the supercurrent, c is the speed of light, and e is the electron charge. In 
our case Vs increases with the distance from the axis of the filament approximately 
proportionally to r. The velocities Vs arising from the neighboring filaments have 
opposite directions in the space between them. Because Vg cannot have discontinuities 
in the superconducting phase, the filaments cannot merge, but must always be divided 
by a boundary where the order parameter "0 = 1^. The situation in the boundary 
region is very similar to that arising in the center of the Abrikosov vortex where 
also changes its sign, demanding that V' = along the vortex axis. 




Figure 3. Profiles of the order parameter for different values of i/'(0) and 
Ho/Hc2 = 0.9. In this case, as well as for a single filament, these profiles are 
independent of ft for 10 < k < 100. The inset shows the total free energy of the 
system of filaments F = nFfu as a function of their density. The solid line is a 
guide to the eye. 



Because the number of filaments in the sample is only limited by a short-range 
repulsion, they are expected to form a rather dense triangular configuration as is 
illustrated in the inset to Fig 4(b). In order to analyze the resulting profile of an 
individual filament, we again use Eqs. (1) and (2) but, in order to satisfy the condition 
tjj = between the filaments, the boundary condition dTp / dr\r^oo = is abandoned. 
Several solutions of the G-L equations are shown in Fig. 3. In this case, the order 
parameter vanishes with a non-zero derivative and we may introduce the radius of the 
filament as the value of r = at the point where — 0. 

As may be seen in Fig. 3, the diameter D of the filament decreases with decreasing 
value of ipiO). The important consequence of the G-L equations is that the diameter 
D of the filament cannot be smaller than a certain minimal value Dmin- The only 
solution of the G-L equations for D < Dmin is = 0. The dependence of Dmin on 
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(1 — Hq/Hc2) is shown in Fig. 4a. We note that also for superconducting lamellae a 
minimal thickness S^ain exists. The dependence of S^ain on the applied magnetic field 
is shown in Fig. 4a, as well. Both -Dmin and (5min decrease with decreasing magnetic 
field, ft has been shown that 6min\Ho=o = t^(,{T) (|]- The corresponding value of 
|i/o=o ~ 4.81^(r). 




1 - HIH,o 1 - HIH, 



Figure 4. (a) The minimal tiiickness of superconducting lamellae, &uih-i) the 
minimal diameter of the superconducting filament, -Dniin, and the equilibrium 
value Deq versus (1 — Hg/Hc'z)- The period d of the vortex structure is shown as 
the solid line. The dashed lines are guides to the eye. The Gibbs free energy F 
of the system of filaments (per unit of volume) with respect to that of the normal 
metal in a magnetic field equal to Hq. The inset represents the triangular lattice 
of superconducting filaments 

In the following we assume that the distance between the filaments is equal to 
their diameter D as is shown in the inset to Fig. 4(b). For a triangular lattice, the 
density n of filaments is thus n = 2/(^/30"^). In order to evaluate the equilibrium 
density of filaments, we plot the total free energy of the system of filaments F — nFju, 
where F is the Gibbs free energy calculated for solutions of the type shown in Fig. 
3, versus n as is illustrated in the inset to Fig. 3. The position of the minimum on 
this curve corresponds to the equilibrium value of n. The magnetic field dependence 
of the equilibrium diameter Deq of the filaments is presented in Fig. 4(a). As may be 
seen, Deq, which is the period of the triangular lattice of the filaments, is considerably 
larger than the equivalent quantity d for the vortex structure. 

The dependencies of tp{0) and (1 — H{0)/Ho) on the magnetic field for an 
equilibrium filament configuration practically coincide with those calculated in the 
single-filament approximation for Hq/Hc2 > 0.5 (see Figs. 2(a) and 2(b)). 

Fig. 4(b) shows the free energy of the system of filaments as a function of 
(1 — Ho/Hc2)- The energy was calculated assuming that the filaments are cylindrical. 
We have to adopt this approach, in order to use one-dimensional G-L equations which 
are the basis of our consideration. It is obvious, however, that, because of their mutual 
interaction, the filaments should adopt a hexagonal rather than a circular cross-section. 
We note, however, that this simplification may only result in free energy values that 
are slightly higher than the actual ones. 

As we have argued, the system of quantized vortices is not the only way to 
realize the mixed state of type-II superconductors which may also be established by an 
ensemble of superconducting filaments. At low magnetic fields, where the equilibrium 



Model of the mixed state of type-II superconductors in high magnetic fields 



6 




■ K= 1.5 - vortices 

• if= 1.5 - filaments 

K= 30 

o K= 100 



o 



•o 



o 



o 



o 



10-1 



10-2 



10-1 
l-HJH, 



10° 



Figure 5. Tlie maximum amplitude of the order parameter for the equilibrium 
configuratioii of filaments and for the vortex lattice. The data for vortices is taken 



from Ref. Q 



diameter of filaments Deq > A(r), the filaments are unstable against the formation of 
a vortex line in the center of the filament. In order to identify the stable arrangement 
in higher fields, the free energies for these two configurations should be compared. 
To our knowledge the only study in which the two-dimensional G-L equations have 
been solved for the vortex lattice is presented in Ref. Q. Fig. 5 shows the results 
of our calculations of the order parameter amplitude V'max versus (1 — Ho/Hc2) for 
filaments, together with the results of Ref. In the high-magnetic-field limit, the 
term proportional to the magnetic moment in the Gibbs free energy may be neglected 
and the density of the free energy reduces to F « —ip'^{r)H'^/8TT. As may clearly be 
seen in Fig. 5, the amplitude of ip for the vortex lattice is considerably smaller than 
that for the filaments and the ratio ?/',jiax(fila'iiients)/'!/'niax(vortices) increases with 
increasing magnetic field. Thus, in sufficiently high magnetic fields, the free energy 
for the configuration of superconducting filaments is expected to be lower than that 
for the vortex lattice. 

The properties of a mixed state consisting of superconducting filaments are quite 
different from those of the Abrikosov vortices. First, because the filaments are 
always separated by normal conducting regions, the sample resistance for currents 
perpendicular to the direction of the magnetic field never vanishes and the true zero- 
resistance superconducting state may be achieved only after the transition to the 
vortex structure. Second, in the case of filaments, the magnetic fiux faces no barriers 
to move in or out of the sample and the magnetization of the sample must be reversible, 
independent of whether the filaments are pinned or not. 

The analysis presented in this paper shows that, in magnetic fields close to Hc2, 
the mixed state may well consist of a triangular lattice of superconducting filaments 
separated by regions where the superconducting order parameter = 0. With 
decreasing external magnetic field, the configuration of superconducting filaments 
necessarily has to undergo a transition to the conventional mixed state, involving 
Abrikosov vortices. The value of the transition field is determined by the free-energy 
balance between these two configurations which cannot be determined without more 
precise calculations of the free energy for both cases. The transition from one type 
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of mixed state to the other involves a complete change of topology and must be 
accompanied by discontinuities in both the resistivity and the magnetic moment of 
the sample. We also expect some hysteresis, as well as a latent heat, dictated by 
the discontinuity of the magnetization. In other words, this transition is expected to 
exhibit all the features of a first-order phase transition. Such transitions are observed 
in high temperature superconductors aX H < Hc2 and are usually attributed to the 
melting of the vortex lattice @, |, |9[ ^ . 

We wish to thank R. Monnier for numerous stimulating discussions. 
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